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1  a  ABSTRACT 200i«»de) 

This  report  describes  the  accomplishments  during  the  contract  period  (March  27,  1990 
to  March  26,  1991)  on  the  computer-aided  analysis  of  CMOS  device  and  circuit  degrada¬ 
tion  due  to  hot-electron  effects.  The  task  involved  four  subtasks:  (1)  modeling  of 
the  gate  oxide  degradation  in  n-channel  MOS  transistor;  (2)  modeling  of  n-channel  MOS 
transistor  behavior  with  localized  oxide  damage;  (3)  simulation  of  gate  oxide  degrada¬ 
tion  during  long-term  circuit  operation;  (4)  determination  of  overall  circuit  perform¬ 
ance  after  hot-electron  stress.  For  the  modeling  of  the  gate  oxide  degradation  in  nMOS 
transistor,  a  localized  triangular  charged  density  distribution  function  has  been  in¬ 
troduced  in  the  drain  end  of  the  channel.  This  model  was  effective  in  explaining  the 
local  electric  potential  near  the  drain,  especially  on  the  flatband  voltage  and  also 
the  changes  in  the  local  channel  electron  mobility.  To  model  the  behavior  of  nMOS 
transistor  with  local  oxide  damages,  a  new  one-dimensional  I-V  equation  has  been  de¬ 
rived  and  implemented  in  a  generic  circuit  simulator,  iSMILE  program.  For  the  modeling 
of  gate  oxide  degradation  process,  the  rate  equations  governing  the  generation  of  in¬ 
terface  states  have  been  implemented  for  dynamic  circuit  operating  conditions.  With 
these  mechanisms  implemented,  the  iSMILE  program  has  been  able  to  simulate  the  degra¬ 
dation  of  circuit  performances  dynamically.  For  large  scale  reliability  simulation, 
sTmnlpr  mnrfels  havp  been _ devi.aed _ and  used  in  the  iDSIM2  program . _ 
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RNAL  REPORT  EVALUATION 
Contract  F30602-88-D-0028  Task  N-0-5701 
VLSI  Design  For  Reliability  -  Hot  Electrons 


The  growth  in  the  use  of  custom,  semi-custom,  and  application  specific  ICs  (ASICs)  in  AF 
systems  has  brought  several  changes  to  the  reliability  qualification  process.  This  occurred 
because  the  costs  associated  with  reliability  characterization  have  become  a  major  part  of 
the  device  cost  when  there  are  only  a  few  hundred  or  thousand  of  such  devices  produced. 
Further,  the  cost  of  a  second  or  third  pass  at  the  design/fabricate/test/evaluate  reliability 
cycle  is  very  time  consuming  as  well  as  expensive. 

Analyzing  an  1C  design  for  susceptibility  to  life  limiting  failure  mechanisms  is  one  approach 
to  enhance  the  probability  that  at  the  end  of  the  first  pass  of  the 
design/fabricate/test/evaluate  reliability  cycle,  the  device  will  have  an  acceptable 
reliability. 

This  effort  addresses  the  development  of  simulation  methods  for  hot  electron  degradation  in 
CMOS  ICs.  It  uses  a  realistic  charge  distribution  in  the  gate  oxide  to  successfully  explain  the 
local  flat  band  voltage,  electric  potential  near  the  drain,  and  the  channel  electron  mobility. 
A  one  dimensional  model  suitable  for  LSI  circuit  analysis  has  been  developed  and  applied  to 
the  simulation  of  a  circuit  with  several  hundred  transistors.  Future  work  will  include 
extending  this  modeling  effort  into  a  macromodeling  scheme  to  address  circuits  with  10K  to 
100K  transistors. 


Other  efforts  being  sponsored  by  Rome  Lab  in  the  area  of  design  analysis  for  ICs  address  the 
estimation  of  current  and  current  density  from  statistical  simulations.  These  estimations 
are  being  applied  to  the  analysis  of  the  susceptibility  of  the  design  to  electromigration  and 
short  term  design  concerns  such  as  IR  drop  and  ground  bounce  in  the  power  and  ground 
busses. 
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VLSI  DESIGN  FOR 
RELIABILITY-HOT  ELECTRON 


Final  Report 

Account  No.  93928-30/Order  No.  RI-68909X 
Contract  F30602-88-D-0028(Task  N-0-5701) 

I.  Introduction 

Long-term  reliability  of  MOS  VLSI  circuits  is  becoming  an  important  issue  as  the  den¬ 
sity  of  VLSI  chips  increases  with  shrinking  design  rules.  In  the  conventional  computer-aided 
design  approach,  physical  aging  of  integrated  circuit  elements  due  to  degradation  mechanisms 
has  been  a  secondary  concern,  and  the  qualification  of  circuit  reliability  has  been  accom¬ 
plished  only  by  accelerated  bum-in  tests  after  manufacture.  This  approach  has  several  draw¬ 
backs  with  respect  to  translating  the  reliability  data  into  physical  design  improvements  for 
long-term  reliability  of  the  circuit.  First,  the  repetition  of  the  "design-manufacture-test"  cycle 
for  incremental  reliability  improvements  is  prohibitively  expensive  in  most  small-volume  pro¬ 
ductions.  Moreover,  some  of  the  long-term  degradation  mechanisms  cannot  be  correctly 
identified  by  accelerated  bum-in  tests.  As  a  result,  the  efforts  to  improve  the  long-term  relia¬ 
bility  of  MOS  VLSI  circuits  have  been  focused  on  process  and  technology  modifications, 
rather  than  design  modifications  of  the  circuit.  The  development  and  use  of  accurate  reliabil¬ 
ity  simulation  tools  is  therefore  cmcial  for  early  assessment  and  improvement  of  circuit  relia¬ 
bility. 

The  hot-carrier  induced  degradation  of  MOS  transistor  characteristics  is  one  of  the  pri¬ 
mary  mechanisms  affecting  the  long-term  reliability  of  MOS  VLSI  circuits.  It  is  likely  to 
become  even  more  critical  in  the  future-generation  chips,  since  the  scaling  of  transistor  dimen¬ 
sions  without  a  corresponding  scaling  of  the  power  supply  voltage  (so  called  constant-voltage 
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scaling)  aggravates  this  problem  [1],[2].  The  hot-carrier  induced  degradation  of  MOS  transis¬ 
tors  is  caused  by  the  injection  of  high-energy  electrons  and  holes  into  the  gate  oxide  region 
near  the  drain.  The  "damage"  is  in  the  form  of  localized  oxide  charge  trapping  and/or  inter¬ 
face  trap  generation,  which  gradually  builds  up  and  permanently  changes  the  oxide-interface 
charge  distribution  [3],[4].  The  extent  of  the  hot-carrier  damage  each  transistor  experiences  is 
determined  by  its  terminal  voltage  waveforms,  i.e.,  by  the  operating  conditions  of  the  circuit. 
Consequently,  the  mechanism  of  hot-carrier  induced  device  degradation  must  be  examined 
within  the  context  of  circuit  simulation. 

The  performance  degradation  in  the  devices  leads  to  the  degradation  of  circuit  perfor¬ 
mance  over  time.  The  reliability  simulation  should  determine  how  the  overall  circuit  perfor¬ 
mance  is  affected  as  a  result  of  device  aging,  and  which  devices  are  most  likely  to  cause  criti¬ 
cal  circuit  performance  failures.  The  framework  of  a  hot-carrier  reliability  simulation  tool 
involves  (i)  modeling  the  gate  oxide  degradation  in  the  MOS  transistor  as  a  function  of  its 
operating  conditions,  (ii)  modeling  the  behavior  of  the  MOS  transistor  with  localized  oxide 
damage,  (iii)  simulation  of  gate  oxide  degradation  during  long-term  circuit  operation,  and  (iv) 
determination  of  the  overall  circuit  performance  after  hot-carrier  stress.  The  integration  of 
these  components  into  a  simulation  framework  is  the  prerequisite  for  accurate  prediction  of 
the  long-term  circuit  reliability  characteristics. 

For  circuit  simulation  purposes,  the  behavior  of  damaged  transistors  was  represented  by 
one-dimensional  homogeneous  models,  based  on  the  assumption  that  the  oxide-interface 
charge  is  uniformly  distributed  over  the  entire  channel  region.  In  their  earlier  work,  Hsu  and 
Tam  have  shown  that  under  the  uniform  charge  profile  assumption  the  amount  of  threshold 
voltage  change  is  linearly  proportional  to  the  amount  of  total  charge  density  [5].  However, 
even  the  definition  of  the  MOSFET  threshold  voltage  loses  much  of  its  physical  meaning  for  a 
channel  with  nonuniform  oxide  charge  and  interface  trap  distribution  [3].  Thus  the  effect  of 
localized  hot-carrier  damage  on  the  I-V  characteristics  of  transistors  cannot  be  accurately 
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accounted  for  by  uniform  device  models. 

The  need  to  represent  the  damaged  MOS  transistor  by  a  nonuniform  and  non- 
symmetrical  structure  has  been  recognized  earlier  [1]  and  several  two-dimensional  modeling 
approaches  have  been  proposed.  Haddara  and  Christoloveanu  have  presented  a  two- 
dimensional  model  for  the  locally  damaged  MOSFET  using  two  channel  regions  and  uniform 
charge  distribution  within  the  damaged  region  [6].  Schwerin  et  al.  have  modeled  the 
behavior  of  tl:e  MOS  transistor  with  fixed  oxide  charges  as  well  as  with  fast  interface  states 
located  near  the  drain  using  the  2-D  device  simulator  MINIMOS  3  [7].  Their  simulation 
results  have  shown  that  the  fast  interface  states  are  necessary  for  accurate  modeling  of  the 
subthreshold  current.  The  nonuniform  fixed  oxide  charge  model  has  also  been  implemented  by 
Roblin  et  al.  for  the  2-D  device  simulator  PISCES  [8]. 

Although  two-dimensional  device  simulation  can  accurately  represent  the  effect  of  local¬ 
ized  hot-carrier  damage  on  the  I-V  characteristics,  it  cannot  be  applied  efficiently  to  simulate 
circuits  with  a  large  number  of  transistors.  For  efficient  circuit-level  simulation,  a  one¬ 
dimensional  device  model  capable  of  representing  the  linear  region  operation  of  the  damaged 
nMOS  transistor  has  been  proposed  earlier  by  the  authors  [9].  In  this  project,  a  new  one¬ 
dimensional  model  is  developd  for  nMOS  transistors  with  hot-carrier  induced  oxide  damage. 
This  MOSFET  model  accounts  for  the  localization  of  the  oxide  damage  near  the  drain  by 
using  a  realistic  charge  density  distribution  profile.  It  is  capable  of  accurately  representing  the 
observed  behavior  of  the  locally  damaged  transistor  both  in  the  linear  and  the  saturation 
regions,  a  task  which  could  not  be  accomplished  using  previous  one-dimensional  models  with 
uniform  oxide  charge  distribution  [10].  The  derivation  of  the  model  is  consistent  with  the 
well-known  gradual  channel  approximation  theory.  The  model  has  been  implemented  in  the 
circuit  simulation  program  iSMILE  [11]  to  enable  efficient  reliability  simulation  of  CMOS  cir¬ 
cuits.  The  model  representation  of  the  hot-carrier  induced  oxide  damage  is  discussed  and  the 
underlying  assumptions  for  the  transistor  model  are  presented  in  Section  II.  The  analytical 
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derivation  of  the  drain  current  model  is  given  in  Section  III,  including  the  modification  of  the 
model  to  account  for  channel  velocity  limitations  in  shon-channel  devices.  The  application  of 
a  pseudo-two-dimensional  approach  to  the  modeling  of  damaged  nMOS  transistors  will  also 
be  discussed.  The  accuracy  of  the  new  model  will  be  demonstrated  in  Section  V  by  compar¬ 
ing  the  simulation  results  with  experimental  data.  The  simulation  of  hot-carrier  induced  dev¬ 
ice  degradation  upon  circuit  performance  are  investigated  in  Section  IV  using  the  iSMILE  cir¬ 
cuit  simulator  and  iDSIM  program.  The  conclusions  will  be  summarized  in  Section  V. 

11.  Hot-Carrier  Induced  Gate  Oxide  Degradation 

A  thorough  understanding  of  the  processes  which  are  responsible  for  the  long-term 
degradation  of  device  characteristics  is  naturally  the  first  prerequisite  for  building  accurate 
simulation  models,  and  for  evaluating  the  device-level  and  circuit-level  damages  associated 
with  these  processes. 

The  hot-carrier  induced  damage  in  nMOS  transistors  is  assumed  to  result  either  in  trap¬ 
ping  of  carriers  on  defect  sites  in  the  oxide  or  the  creation  of  interface  states  at  the  silicon- 
oxide  interface,  or  both.  Recent  studies  have  shown  that  the  charge-pumping  technique  can 
provide  precise  and  direct  information  on  the  interface  properties  in  MOS  transistors  [3], [12]. 
In  addition,  this  technique  is  capable  of  independently  providing  the  amount  of  interface  traps 
generated  during  injection,  and  the  amount  of  charges  that  have  been  trapped  in  the  gate 
oxide,  even  for  the  case  of  localized  degradation. 

Two  distinct  voltage  regions  of  stressing  have  been  identified  which  indicate  different 
types  of  degradation.  For  the  region  of  maximum  substrate  current  (V^s  =  2  Vgs).  the  oxide 
damage  has  been  linked  to  interface  trap  generation  through  hot-electrons  and  hot-holes,  and 
no  significant  charge  trapping  could  be  observed  [2]-[4].  For  the  region  of  high  electron 
injection  into  the  oxide  (Vos  =  Ves),  the  damage  is  caused  to  a  great  extent  by  charge  trap¬ 
ping,  and  also  by  moderate  interface  trap  generation  [41,[8].  Consequently,  an  investigation  of 
hot  carrier  induced  degradation  mechanisms  must  encompass  the  charge  trapping  and  the 
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interface  trap  generation  in  nMOS  transistors,  both  of  which  are  caused  by  hot-carriers 
injected  into  the  gate  oxide. 

A.  Charge  Trapping 


When  the  transistor  is  operating  in  or  near  the  saturation  region,  the  electrons  moving 
from  the  source  to  the  drain  experience  a  significant  velocity  increase  due  to  the  high  horizon¬ 
tal  electric  field.  A  small  percentage  of  these  high-energy  electrons  are  deflected  towards  the 
gate  electrode  by  elastic  scattering.  If  the  kinetic  energy  of  the  deflected  electron  (hot  elec¬ 
tron)  is  sufficient  to  overcome  or  tunnel  through  the  potential  energy  barrier  at  the  silicon- 
oxide  interface,  the  electron  is  injected  into  the  gate  oxide  near  the  drain.  This  process  can  be 
characterized  by  the  electron  injection  current  into  the  gate  oxide  in  terms  of  drain  current  Iq 
and  other  parameters; 


*EI 
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is* 

2 

K 

. ‘J>b  . 

P(Eox)  exp 

(1) 


Here,  denotes  the  maximum  channel  electric  field  in  the  direction  of  the  channel  current, 
<l>b  denotes  the  Si-Si02  potential  barrier  height  and  P(Eox)  denotes  the  probability  that  an  elec¬ 
tron  can  overcome  this  potential  barrier.  It  is  assumed  that  this  injection  is  confined  to  a  rela¬ 
tively  small  portion  of  the  channel  near  the  drain.  The  electron  injection  current  reaches  its 
peak  value  for  Vqs  =  Vps.  which  is  identified  as  a  bias  condition  favorable  for  charge  trapping 
14],  The  relationship  between  the  injected  electron  current  and  the  trapped  oxide  charge  den¬ 
sity  Nox  can  be  described  by 


No*(Ni„j)  =  N,  [l  -  e'  +  Nj  [l  -  e"  .  (2) 

where  Oj  and  02  represent  the  trapping  cross-sections  and  the  injected  electron  fluence  Ni„j  is 
calculated  as  the  integral  of  electron  injection  current  density  over  time  [9]. 
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B.  Interface  Trap  Generation 

The  other  prominent  hot-carrier  injection  mechanism  is  based  on  impact  ionization,  in 
which  high-energy  electrons  moving  horizontally  along  the  channel  create  elecnon/hole  pairs 
by  collisions  near  the  drain.  These  hot  carriers  have  been  identified  as  being  responsible  for 
the  generation  of  oxide-interface  states  (traps),  a  process  which  is  distinctly  different  from 
charge  trapping.  The  hot-electron  induced  trap  generation  process  can  be  characterized  by  the 
bond-breaking  current  Igg  [2]: 

Ibb  =  C,  loexp  — (3) 
.Q  ^  . 

where  represents  the  critical  energy  for  elections  to  create  interface  traps,  and  is  the 
mean  free  path  for  electrons.  A  similar  expression  can  be  derived  for  the  hot-hole  component 
of  Ibb.  Using  the  rate  equations  governing  the  generation  of  interface  states  (N^),  the  long¬ 
term  interface  charge  accumulation  under  dynamic  operating  conditions  is  described  by  the 
following  simple  differential  equation  [20]. 

dN:,  r  Xh 

— i  l+Bp~Nii  =K<Ibb>  (4) 

dt  L  Dh  J 

Here,  <Ibb>  represents  the  average  value  of  the  bond-breaking  current  over  a  period,  Dh  and 
Xh  represent  diffusion  constants,  K  and  Bp  represent  process-dependent  coefficients. 

The  charge  trapping  model  given  by  Eqs.  (1)  and  (2),  and  the  interface  trap  generation 
model  given  by  Eqs.  (3)  and  (4)  can  be  applied  for  predicting  the  long-term  dynamics  of  hot- 
cairier  induced  device  degradation.  It  must  be  noted  that  the  effects  of  hot-carrier  degradation 
upon  current-voltage  characteristics  of  nMOS  transistors  cannot  be  explained  by  simple 
parameter  changes  in  the  conventional  MOSFET  simulation  models.  To  account  for  the 
behavior  of  the  damaged  nMOS  transistors  using  a  conventional  SPICE  model  such  as  the 
BSIM,  most,  if  not  all,  of  the  model  parameters  must  be  modified  accordingly  [21].  Even 
such  an  approach  would  be  largely  empirical,  since  the  conventional  uniform  MOSFET 
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models  do  not  accommodate  localized  modifications  of  physical  device  parameters.  This  fact 
complicates  the  simulation  of  the  damaged  device  behavior  with  various  levels  of  hot-carrier 
damage. 

The  localization  of  hot-carrier  injection  into  the  gate  oxide  leads  to  the  localization  of  the 
resulting  oxide  damage  in  the  form  of  charge  trapping  and  interface  trap  generation 
[3], [4], [22].  The  change  of  oxide  characteristics  due  to  charge  injection  occurs  in  a  relatively 
short  channel  section  (approximate  length  of  0.1  pm)  at  the  drain  end.  This  localization  of 
damage  largely  complicates  the  interpretation  of  damaged  transistor  characteristics,  and  it  is 
also  responsible  for  the  observed  asymmetry  of  current-voltage  characteristics  in  damaged 
nMOS  transistors. 

Although  it  is  quite  clear  that  the  hot-carrier  induced  oxide  charge  is  localized  very  near 
to  the  drain,  the  shape  of  the  charge  distribution  profile  is  not  known  exactly.  The  existing 
measurement  results  indicate  that  the  peak  value  of  the  charge  profile  and  the  length  of  the 
damaged  region  increase  with  hot-carrier  stress.  Haddara  et  al.  [6]  and  Roblin  et  al.  [8]  have 
used  uniform  (constant)  charge  profiles  in  two-dimensional  device  simulations,  whereas 
Schwerin  et  al.  have  assumed  a  Gaussian  charge  distribution  [7].  In  all  cases  the  length  of  the 
damaged  region  is  assumed  to  be  about  0.1  pm.  Ancona  et  al.  [12]  have  reported  charge 
pumping  measurement  results  which  indicate  a  triangular  charge  distribution  profile  with  a 
sharp  peak  very  close  to  the  drain  (Fig.  1).  For  the  derivation  of  model  equations,  the  local¬ 
ized  oxide-interface  charge  due  to  hot-carrier  injection  will  be  represented  by  a  simple  tri¬ 
angular  charge  density  distribution  profile  (Fig.  2). 


Here,  Li  denotes  the  length  of  the  undamaged  channel  region,  whereas  L2  denotes  the  length 


of  the  damaged  region.  Notice  that  the  amount  of  oxide-interface  damage  is  described  by 
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only  two  parameters,  i.e.,  Qpe*k  and  L2.  The  hot-carrier  induced  oxide-interface  charge  has  a 
significant  influence  on  two  device  properties:  (i)  it  affects  the  local  electric  potential  in  the 
adjacent  drain  region,  and  (ii)  it  changes  the  local  channel  electron  mobility.  The  potential 
effect  can  be  accounted  for  as  an  additional  term  in  the  flat-band  voltage  (Vfb)  expression. 


VpaCy)  =  **\fs  "* 


Qii(y) 


(6) 


Here  <I>ms  represents  the  work  function  difference,  Q,,,  represents  the  constant  positive  oxide- 
interface  charge  density  in  C/cm^  which  is  independent  of  hot-carrier  i.ijection,  and  Cq, 
represents  the  gate  oxide  in  F/cm^.  The  reduction  of  channel  electron  mobility  can  be 
represented  by  Mathiessen’s  rule  as  proposed  by  Nishida  and  Sah  [13].  This  approach 
assumes  that  the  total  reciprocal  mobility  is  composed  of  bulk  mobility  Po  ^^d  oxide-interface 
charge  scattering  mobility  Po,  as  follows. 


m,(y)  m)  My) 

The  electron  mobility  due  to  charge  scattering  is  given  by 

Pox(y)  =  1000  (3xl0“/Qi,(y))  (T/80)  [cm^/V  s]  (8) 

where  Qii(y)  is  the  area-density  of  the  oxide-interface  charge.  The  derivation  of  drain  current 
model  equations  is  presented  in  the  following  section. 


III.  Modeling  of  nMOS  Transistor  with  Gate  Oxide 
Damages  for  Reliability  Simulation 

The  drain  current  model  derivation  is  based  on  the  assumption  that  the  gradual-channel 
approximation  is  valid  for  the  damaged  nMOS  transistor.  This  approximation  states  that  the 
electric  field  in  the  direction  of  current  flow  is  much  smaller  than  the  field  perpendicular  to 
the  silicon  surface,  thus  allowing  the  one-dimensional  analysis  of  the  drain  current  Iq. 
Throughout  this  section,  all  voltage  values  will  be  given  with  respect  to  the  source  potential, 
unless  stated  otherwise.  The  mobile  electron  charge  Qn  along  the  channel  region  is  expressed 
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by 


Q„(y)  =  -  Co,  [vG-Vc(y)-VFB(y)-2l<I>pl  ]  +  [2esqN.  (2l<I>pl+Vc(y)-VB)]  ^  (9) 

where  Vq,  Vb  and  Vc  represent  the  gate  voltage,  the  substrate  voltage  and  the  channel  voltage, 
respectively.  First,  the  derivation  will  be  carried  out  in  the  linear  operation  region,  where  a 
continuous  channel  is  assumed  to  exist  between  the  source  and  the  drain. 


A.  Linear  Region  Equations 

When  the  nMOS  transistor  is  operating  in  the  linear  region,  an  incremental  length,  dy, 
along  the  channel  sustains  a  voltage  drop,  dVc,  that  can  be  expressed  as  the  well-known  pro¬ 
duct  of  the  drain  current  Iq  and  the  incremental  channel  resistance; 


or 


dVc(y)  =  -  Id 


_ dy _ 

W  p„(y)  Q„(Vc.y) 


(10) 


Id  dy  =  -  W  p„(y)  Q„(y)  dVc  .  (11) 

The  solution  of  this  differential  equation  will  be  substantially  simplified  using  the  following 
observations.  First,  the  oxide-interface  charge  distribution  profile  Qit(y)  is  equal  to  zero 
between  y=0  and  y=Li  according  to  Eq.  (5).  Thus  the  right-hand  side  of  Eq.  (11)  is  only  a 
function  of  the  channel  voltage  in  this  undamaged  region.  To  achieve  a  similar 
simplification  in  the  damaged  region  of  length  L2,  the  charge  density  will  be  assumed  to  be 
constant  and  equal  to  the  average  value  of  Qi,(y)  in  this  region,  i.e.,  Qu=Qpeak/2.  This  assump¬ 
tion  actually  reduces  the  oxide-interface  charge  distribution  given  in  Fig.  2  to  a  locally  uni¬ 
form  profile  for  the  linear  operation  region.  The  benefit  of  using  the  triangular  charge  distri¬ 
bution  profile  will  be  demonstrated  for  the  saturation  region  in  sub-section  B.  A  uniform 
charge  distribution  profile  in  the  damaged  region  implies  that  the  channel  electron  mobility  in 
this  region  (p,^)  is  independent  of  y  and  that  Q„  is  only  a  function  of  Vq.  Under  these  simpli¬ 
fying  assumptions,  the  solution  of  Eq.  (11)  can  be  found  in  two  adjacent  channel  regions. 
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The  parameters  associated  with  the  undamaged  region  (Region  1)  and  the  damaged  region 
(Region  2)  will  be  indicated  below  with  the  subscript  indices  1  and  2,  respectively.  Integrat¬ 
ing  both  sides  of  Eq.  (11)  from  the  source  to  the  boundary  of  the  damaged  region  (Li)  yields 


Id  Jdy  =  -W  J  My)  Qn(y)  dVc 


=  Wm.,Co»|  Vg-Vpb,-2I<1»pI-Vc-^  [2esqN.(2IC»pl+Vcr-VB)]^  dVc . 

0  I 


Id  =  —MCo,  -  (Vg-Vfb,-2I<I>pI)  Vp  - 

Lj  I  I 


2  V2£sqN,  4  4 

f  ^  *  (2l<I>pl-VB+Vp)2 -(2I<I>pI-Vb)2 

J  w/w  L 


In  Eqs.  (12)  and  (13),  Vp  represents  the  channel  voltage  at  the  boundary  of  the  damaged 
region,  i.e.,  at  y=Li.  It  is  easily  seen  that  Eq.  (13)  corresponds  to  the  exact  drain  current 
equation  of  an  undamaged  transistor  with  channel  length  Lj  and  drain  voltage  Vp.  Since  the 
value  of  Vp  is  not  known  yet,  Eq.  (II)  must  now  be  solved  in  the  damaged  region  with  the 
corresponding  boundary  conditions.  The  integration  from  the  boundary  of  the  damaged  region 
(y=Li)  to  the  drain  is  carried  out  as  follows. 


Id  J  dy  =  -W  I  My)  Qn(y)  dVc 


=  Wp„2CoxJ  Vg-Vfb2-'2I<I>pI-Vc-7^  [2esqN,(2l<I)pl-^Vc-VB)l^  dV^ 
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w 

Id  =  7“mi2Cai 


*  (Vo-VpBr-Z  1  <*>?  I )  (Vo-Vp)  -  j-  (Vg-V|) 


2  V2£sqN. 

3 


1  1 
(2  l<I>p  I  -Vb+Vd)  *  -  (2 1  Op  I  -Vb+Vp)  2 


(15) 


With  some  simple  manipulation,  Eq.(lS)  can  be  written  in  the  following  form. 

(Vd-Vp)2 


Id  =  ■^mi2Co»  1 


[(Vcr-Vp)-VpB2-2l<I)pl  ]  (Vo-Vp)  - 


2  V2esqN, 

3 


[2IOpl-(VB-Vp>f(VD-Vp)]^  -  jilOpKVB-Vp)] 


(16) 


Upon  inspection,  Eq.  (16)  can  be  identified  as  the  exact  drain  current  equation  of  an 
nMOS  transistor  with  channel  length  L2,  drain  bias  Vq  and  source  bias  Vp.  Thus,  the  opera¬ 
tion  of  the  damaged  nMOS  transistor  in  the  linear  region  is  represented  by  two  transistors 
connected  in  series.  The  drain  current  value  Iq  for  given  Wq,  Vp  and  Vb  can  be  obtained  by 
solving  Eqs.  (13)  and  (15)  with  Newton-Raphson  iteration.  The  onset  of  conduction  in  the 
channel  is  also  determined  iteratively,  since  the  threshold  voltage  of  the  locally  damaged 
nMOS  transistor  cannot  be  represented  by  a  closed-form  expression.  However,  the  following 
empirical  approximation  for  the  threshold  voltage  has  been  found  to  be  sufficiently  accurate. 

V-n.  =  Vfb2  +  2l<I»pl  +  ^  ^  V2T^^i^.  (17) 

Here,  VpB2  represents  the  average  flat-band  voltage  of  the  damaged  channel  region. 

B.  Saturation  Region  Equations 

The  analysis  leading  to  the  current  equations  (13)  and  (IS)  is  valid  as  long  as  the  inver¬ 
sion  charge  exists  along  the  entire  length  of  the  channel.  When  the  drain  voltage  Vq  is 
increased  sufficiently  to  cause  the  mobile  electron  charge  Q„(L)  at  the  drain  to  decrease  to 
zero,  the  channel  pinches  off  and  the  drain  current  saturates.  Notice  that  this  definition  of 
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current  saturation  does  not  take  into  account  the  channel  velocity  saturation  effect  which  dic¬ 
tates  the  saturation  conditions  in  short  channel  transistors,  and  should  therefore  be  viewed  as  a 
first-order  approximation.  The  modifications  according  to  channel  velocity  limitations  will  be 
covered  in  the  next  section. 

The  drain  voltage  corresponding  to  the  onset  of  saturation  can  be  obtained  by  setting 
Q„(L)  =  0  in  Eq.(9)  and  solving  for  Vc(L)  =  Vdsat- 


'^DSAT  =  Vq-V  fb(L>-2  I  <I>p 


^  esqN,  J 


1+- 


2CL 


EsqN. 


(Vg-Vpb(L>-Vb)  h 


(18) 


Note  that  the  flat-band  voltage  Vpg  in  Eq.(18)  is  calculated  using  the  triangular  oxide  charge 
distribution  profile,  and  not  the  average  value  of  the  oxide  charge  in  the  damaged  region.  As 
the  drain  voltage  Vq  is  increased  beyond  while  the  gate  voltage  is  held  constant,  the 

effective  channel  length,  i.e.,  the  portion  of  the  channel  along  which  surface  inversion  exists, 
decreases  due  to  channel  length  modulation.  In  the  damaged  transistor  model  presented  here, 
the  channel-end  voltage  Vce  beyond  saturation  is  calculated  by  taking  into  account  the  non- 
uniform  oxide  charge  distribution  given  in  Eq.(5).  Thus,  contrary  to  the  conventional  model 
of  the  undamaged  nMOS  transistor  in  saturation,  the  channel-end  voltage  does  not  remain 
equal  to  V^s^t  but  increases  as  the  pinch-off  point  moves  toward  the  source.  For  Vp  >  Vdsat* 
the  saturation  drain  current  of  the  transistor  is  found  by  solving  Eqs.  (13)  and  (15),  where  the 
drain  voltage  Vq  in  the  equations  is  replaced  by  the  channel-end  voltage 


The  length  of  the  depletion  region  between  the  drain  and  the  channel-end  point  is  calcu¬ 
lated  by  combining  the  lateral  electric  field  components  due  to  (i)  the  potential  difference 
between  the  drain  and  channel-end  point,  (ii)  the  fringing  electric  field  between  the  gate  and 
the  drain,  and  (iii)  the  fringing  electric  field  between  the  gate  and  the  channel-end  point.  An 
approximate  formula  for  the  length  AL  of  channel-end  depletion  region  has  been  given  by 
Frohman-Bentchkowsky  and  Grove  [14].  The  formula  is  modified  according  to  the  non- 
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uniform  oxide  charge  distribution  as  follows. 

_L-  /  "  q  N, - ^  ep  g  [Vp-Vc^]  +  p  [Vg(AL)-Vce(AL)] 

AL  “  \  2es  [Vd-Vce(AL)]  EsW  Vq  -  Vce(AL) 

where 


(19) 


Vq(AL)  =  Vc--~  [Qu(L-AL)  -  Qo,]  (21) 

Vg'  =  Vg-“  [Qp«Jc-Qox]  .  (22) 

The  coefficients  a  and  p  represent  fringing  field  factors  accounting  for  the  lateral  contributions 
of  the  fringing  electric  fields.  The  channel-end  depletion  region  length  AL  for  given  bias  con¬ 
ditions  is  calculated  from  Eqs.  (19)-(22)  by  Newton-Raphson  iteration.  Thus  the  nMOS 
transistor  operating  in  saturation  region  can  be  accurately  modeled  by  a  transistor  with 
effective  channel  length  L-AL  and  with  the  corresponding  drain  voltage  Vce(AL).  Linear 
region  current  equations  (13)  and  (15)  are  used  to  compute  the  drain  current  since  this 
modified  transistor  operates  at  the  saturation  boundary.  Note  that  the  average  oxide  charge 
density  used  in  Eqs.  (13)  and  (15)  decreases  as  the  pinch-off  point  moves  away  from  the 
drain,  because  the  oxide  charge  profile  directly  above  the  channel-end  depletion  region  has  no 
effect  on  the  drain  current.  When  the  length  AL  of  the  depletion  region  exceeds  the  length 
of  the  damaged  oxide  region,  the  drain  current  of  the  damaged  transistor  becomes  approxi¬ 
mately  equal  to  the  drain  current  of  an  undamaged  transistor  in  saturation. 

C.  Consideration  of  Channel  Velocity  Limitations 

The  mobile  electron  charge  Q„  moves  in  the  channel  under  the  influence  of  Ey,  the  elec¬ 
tric  field  component  tangential  to  the  silicon-oxide  interface.  At  the  same  time,  Q„  is  held 
near  the  interface  by  the  field  component  E,,  which  is  perpendicular  to  it.  Both  field 
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components  E,  and  Ey  influence  the  velocity  of  the  electrons  moving  along  the  channel. 

Experiments  have  shown  that  carrier  velocities  near  the  surface  tend  to  saturate  at  higher 
tangential  fields  [15],[23].  This  behavior  can  be  modeled  by  the  following  expression 


V  =: 


»m.Eyi 


(1  +  IEy/Ecl®)‘'“ 


(23) 


where  v  represents  the  carrier  velocity,  and  a  and  Ec  are  empirical  parameters.  This  expression 
can  be  further  simplified  by  assuming  a=l  [16].  Using  Eq.  (23)  with  a=l  and 
Ey  =  -  dVc/dy,  Eq.  (11)  can  be  rewritten  as 


m,EyWQ„(y) 

‘D  - - 

1  +(Ey/Ee) 

This  equation  is  reorganized  and  integrated  along  the  channel  to  yield; 

Vd 


(24) 


Id  = - - -  JlinQndVc  .  (25) 

L(1  +  IVd/EcD  i 

Thus  the  influence  of  channel  velocity  saturation  can  be  interpreted  as  a  lengthening  of  the 
channel  by  IVu/EcI,  or  as  a  reduction  of  the  electron  mobility  while  the  channel  length  is 
kept  constant. 


To  approximate  the  drain  voltage  Vqsat  at  the  onset  of  saturation,  the  following  semi- 
empirical  expression  has  been  used. 


VDSAT<new)  =  7  lE^LI 


2  VpSAT 
lEcLI 


(26) 


Here,  the  coefficient  y  with  a  value  between  0.9  and  1.0  accounts  for  the  fact  that  the  channel 
actually  saturates  brfore  the  mobile  electron  charge  Q„  is  reduced  to  zero  at  the  drain. 


In  addition  to  the  effects  of  channel  velocity  saturation,  the  influence  of  the  vertical  elec¬ 
tric  field  E,  upon  carrier  transport  must  also  be  considered.  This  field  component  influences 
the  scattering  at  the  silicon-oxide  interface,  and  causes  a  virtual  reduction  of  the  carrier  mobil¬ 
ity  in  the  channel  [15].  This  effect,  combined  with  the  carrier  velocity  saturation  examined 
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above,  has  been  modeled  empirically  as 

Hrff  - - - -  •  (27) 

1  +  IVc/EcLI  +tiVg 

In  this  expression,  E,.  represents  the  characteristic  field  for  the  channel  velocity,  and  r| 
represents  the  empirical  parameter  corresponding  to  the  vertical  field  effect.  Equation  (27) 
can  be  used  for  both  regions  of  the  damaged  MOSFET;  by  setting  the  drain  voltage  equal  to 
Vp  in  the  undamaged  Region  I,  and  Vq-Vp  in  the  damaged  Region  II.  With  this  simple 
improvement,  the  proposed  MOSFET  model  accurately  simulates  the  behavior  of  transistors 
with  channel  lengths  down  to  1  pm. 

D.  Pseudo  Two-Dimensional  Modeling  of  Damaged  MOSFETs 

The  one-dimensional  model  equations  in  this  section  were  derived  using  the  assumption 
that  the  gradual  channel  approximation  is  valid  for  the  damaged  nMOS  transistor,  and  in  par¬ 
ticular,  for  the  short  damaged  region  at  the  drain  end  of  the  channel.  However,  this  assump¬ 
tion  becomes  questionable  for  transistors  with  very  short  channels,  i.e.,  for  sub-micron 
geometries,  where  the  potential  distribution  becomes  two-dimensional.  Thus,  the  calculation 
of  the  depletion  region  length  and  the  saturation  voltage  requires  the  accurate  estimation  of 
the  vertical  as  well  as  the  horizontal  electric  field  components  in  the  drain  region  using  a 
pseudo  two-dimensional  modeling  approach  [17],[18]. 

When  the  nMOS  transistor  is  operating  in  saturation,  the  channel  is  divided  into  two 
regions:  The  source  region  where  the  Gradual  Channel  Approximation  (GCA)  is  valid,  extend¬ 
ing  from  the  source  to  the  point  of  saturation,  and  the  drain  region  extending  firom  the  point 
of  saturation  to  the  drain  junction.  The  saturation  current  can  be  found  by  setting  the 
channel-end  voltage  equal  to  Vq^at  ^nd  channel  length  equal  to  L-AL  in  the  linear  region 
equations,  and  solving  for  Iq.  Thus,  the  drain  depletion  region  length  AL  and  the  saturation 
voltage  Vqsat  ^ust  be  determined  using  the  pseudo  two-dimensional  analysis.  In  the  drain 
region,  mobile  electrons  are  assumed  to  spread  over  an  average  depth  X^.  A  rectangular 


17 


Gaussian  surface  that  encloses  both  the  mobile  and  bulk  charges  is  shown  in  Fig.  3.  Applying 
Gauss’s  law  to  the  sides  of  the  rectangle  yields  the  following  expression. 

X*,  y  X^ 

-  J  e,Eidx  +  Je,  E2dy+  J  e^Esdx 
0  0  0 

y  y 

-Jeo»E4dy  =  -qJ 
0  0 

Here,  Ei,  Ej,  E3  and  E4  are  the  electric  fields  perpendicular  to  the  boundaries  of  the  Gaussian 
surface.  Note  that  the  field  component  E4,  which  represents  the  oxide  electric  field  above  the 
drain  depletion  region,  also  contains  the  contribution  of  the  localized  oxide  charge  due  to 
hot-carrier  damage.  An  iq)proximate  solution  of  Eq.  (28)  can  provide  the  saturation  voltage 
Vqsat  as  well  as  the  depletion  region  length  AL  as  a  function  of  the  external  bias  conditions 
[18],[19].  However,  this  approach  still  requires  a  simultaneous  solution  of  five  non-linear 
equations,  which  significantly  increases  the  computational  complexity  of  the  model  and  dimin¬ 
ishes  its  applicability  for  circuit-level  simulation  tasks. 

E.  Simulation  Results 

The  accuracy  and  various  properties  of  the  damaged  MOSFET  current  model  will  be 
examined  by  comparing  the  simulation  results  with  experimental  data.  Circuit  simulation 
examples  will  also  be  presented  to  indicate  the  implementation  of  the  new  model  for  simulat¬ 
ing  hot-carrier  induced  circuit  performance  degradation. 

The  measured  and  simulated  Iq-Vq  characteristics  of  an  nMOS  transistor  before  and  after 
hot-carrier  stress  are  given  in  Fig.  4.  The  experimental  data  on  the  damaged  nMOS  transistor 
were  obtained  from  Texas  Instruments,  Inc.  The  geometry  of  the  measured  nMOS  transistor 
is  described  by  the  following  parameters:  W  =  10  pm,  L  =  1.0  pm,  tg,  =  20  ran.  The  transistor 
has  been  stressed  under  maximum  hot-carrier  injection  conditions  (Vp  =  6.5  V,  Vq  =  2.6  V)  for 
423  hours.  The  simulated  Iq-Vd  characteristics  correspond  to  a  peak  interface  charge  density 
of  Qp^  =5  q-8.2xI0"  C/cm^  and  to  a  damaged  region  length  of  L2  =  0.12pm.  The  model 


J  (n  +  NA)dx 


(28) 


dy 


IS 


parameters  which  had  not  been  specified  by  the  device  manufacturer  were  optimized  using  the 
parameter  extraction  routine  in  iSMILE.  A  list  of  model  parameters  and  their  values  are 
given  in  Table  1. 

The  comparison  with  the  measurement  results  demonstrates  that  the  proposed  model  is 
capable  of  accurately  representing  the  behavior  of  the  nMOS  transistor  after  localized  hot- 
cairier  damage.  Figure  5  shows  the  measured  and  the  simulated  output  conductance  of  the 
undamaged  transistor.  For  accurate  circuit  reliability  simulation,  the  transistor  model  must  be 
able  to  represent  the  behavior  of  the  damaged  device  for  different  hot-carrier  damage  levels, 
which  correspond  to  different  stress  times  under  DC  conditions.  In  Figure  6,  the  measured 
drain  current  degradation  versus  hot-carrier  stress  time  is  compared  with  simulation  results. 
The  drain  current  degradation  is  defined  as  the  decrease  of  the  current  at  Vq  =  5  V, 
Vq  =  1.5  V.  The  DC  stress  conditions  are  Vq  =  6.5  V,  Vq  =  2.6  V. 

The  forward  and  reverse  Iq-Vq  characteristics  of  a  damaged  nMOS  transistor  are  given 
in  Fig.  7  [7].  In  this  example,  the  channel  length  of  the  transistor  is  L  =  1.7  pm,  and  the  stress 
conditions  are  Vq  =  8  V,  Vq  =  3  V,  =  5  x  10^  s.  The  proposed  MOSFET  model  replicates 
the  strong  asymmetry  between  forward  and  reverse  operation,  which  is  due  to  the  localization 
of  oxide-interface  damage  near  the  drain.  The  uniform  charge  distribution  model  [2],[5]  is  not 
capable  of  reproducing  this  asymmetry  which  may  seriously  affect  the  performance  of  bidirec¬ 
tional  devices  such  as  transmission  gates.  On  the  other  hand,  the  uniform  model  was  found  to 
provide  satisfactory  simulation  results  for  the  case  where  both  the  source  and  drain  ends  of 
the  channel  are  subjected  to  equal  amounts  of  hot-carrier  damage.  Figure  8  shows  the  simu¬ 
lated  Id-Vg  characteristics  of  an  nMOS  transistor  with  different  levels  of  hot-carrier  damage. 
As  expected,  the  threshold  voltage  increases  with  the  amount  of  oxide-interface  charge. 
Notice  that  the  maximum  channel  transconductance  also  increases  as  a  result  of  localized 
degradation.  This  phenomenon,  which  was  also  observed  in  two-dimensional  device  simula¬ 
tions  with  PISCES  using  fixed  negative  oxide  charges  in  the  drain  region  [8],  is  due  to  the  ini- 
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tial  model  assumption  that  all  generated  interface  traps  are  occupied  with  electrons. 

The  circuit-level  implementation  of  the  proposed  MOSFET  model  using  the  iSMILE  cir¬ 
cuit  simulator  is  demonstrated  in  the  next  example.  A  7-stage  CMOS  ring  oscillator  with  0.2 
pF  loading  at  each  stage  is  considered.  The  (W/L)-ratios  of  the  nMOS  and  pMOS  transistors 
are  (10/2)  and  (15/2),  respectively.  To  investigate  the  influence  of  hot-carrier  induced  device 
damage  upon  the  circuit  performance,  a  localized  oxide-interface  charge 
(Qpeak  *  q-1.5xl0‘^  C/cm^,  L2  =0.1  pm  )  is  introduced  in  all  nMOS  transistors,  and  the  circuit  is 
simulated  with  iSMILE.  Figure  9  shows  a  degradation  (decrease)  of  the  ring  oscillator  fre¬ 
quency  by  approximately  10%  due  to  the  localized  oxide  damage  in  nMOS  transistors.  To 
compare  this  result  with  uniform-damage  model  results,  the  same  amount  of  interface  charge 
is  assumed  to  be  distributed  uniformly  along  the  channel  region  of  each  transistor,  and  the  cir¬ 
cuit  is  simulated  with  iSMILE.  It  is  found  that  the  observed  degradation  of  the  oscillation 
frequency  is  about  5%,  indicating  that  a  transistor  model  based  on  uniform  charge  distribution 
tends  to  underestimate  the  circuit  performance  degradation.  The  simulated  operation  of  a 
cross-coupled  latch  circuit  is  shown  in  Fig.  10  with  (i)  locally  damaged  transistors,  and  (ii) 
uniformly  damaged  transistors.  As  in  the  previous  example,  it  can  be  seen  that  the  uniform- 
damage  assumption  again  underestimates  the  circuit  performance  degradation.  The  rise  time 
and  the  fall  time  of  the  respective  terminal  voltages  after  degradation  are  found  to  be 
significantly  larger  using  the  localized  damage  model. 

Finally,  the  hot-carrier  induced  degradation  of  an  nMOS  transistor  in  bidirectional  opera¬ 
tion  is  examined  in  the  following  example.  Figure  11  shows  a  simple  nMOS  pass-transistor 
structure  in  which  the  capacitances  C]  and  C2  are  charged  and  discharged  during  each  cycle. 
The  pass-transistor  M3  experiences  significantly  more  hot-carrier  degradation  during  the 
charge  phase  than  during  the  discharge  phase  because  the  initial  drain-to-source  voltage  is 
higher  in  the  former  case.  Consequently,  transistor  M3  is  expected  to  sustain  a  unidirectional 
oxide  damage  despite  of  its  bidirectional  operation  mode.  The  result  of  this  one-sided  local 
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damage  is  that  the  drain  current  is  significantly  reduced  during  the  discharge  phase,  since  the 
oxide  damage  is  then  located  at  the  source  end  of  the  channel.  The  simulated  input  and  out¬ 
put  voltage  waveforms  are  shown  in  Figure  12.  It  is  seen  that  the  fall  time  of  the  output  vol¬ 
tage  waveform  increases  with  hot-carrier  induced  damage,  because  the  drain  current  of  M3 
during  discharge  phase  is  substantially  reduced. 

rv.  Reliability  Simulation  of  Very  Large  CMOS  Circuits 

In  the  previous  section,  we  have  presented  a  new  approach  for  simulating  hot-canier 
induced  MOS  circuit  degradation.  However,  the  detailed  circuit  simulation  needed  for  deter¬ 
mining  the  stress  conditions  of  individual  devices  restricts  the  computational  efficiency  of  reli¬ 
ability  simulation  approaches  for  very-large  scale  integrated  circuits. 

Figure  13  shows  a  hierarchical  simulation  approach  presented  for  estimating  the  hot- 
carrier  induced  degradation  of  device  characteristics  and  circuit  performance  in  large  scale 
MOS  circuits.  The  reliability  simulation  tool  presented  here  incorporates  a  two-tier  hierarchi¬ 
cal  simulation  approach  for  improving  the  computational  efficiency.  Large  circuits  are  simu¬ 
lated  using  a  fast  timing  simulator  to  detect  subcircuits  likely  to  cause  reliability  problems. 
Then,  detailed  circuit  simulation  is  performed  on  the  suspected  subcircuits.  An  accurate  one¬ 
dimensional  MOSFET  model  and  a  repetitive  simulation  scheme  are  adopted  to  ensure  accu¬ 
rate  prediction  of  the  circuit-level  degradation  process  under  dynamic  operating  conditions. 
The  fast  simuladon  is  performed  using  the  mixed-mode  simulator  iDSIM2,  whereas  iSMILE 
is  used  for  the  detailed  simulation  [24]. 

iDSIM2  reads  in  a  SPICE-type  circuit  description  and  creates  a  circuit  database.  Based 
on  this  database,  iDSIM2  performs  the  following  preprocessing.  First,  it  partitions  the  MOS 
circuit  into  channel-connected  subcircuits  and  adds  the  information  back  to  the  database. 
Each  subcircuit  may  contain  transistors,  resistors,  diodes  and  capacitors.  It  detects  feedback 
loops  existing  in  the  circuit  based  on  the  basic  principle  that  signals  propagate  from  gate  node 
to  drain  and  source  nodes  of  a  MOS  transistor  while  the  signal  propagation  between  drain  and 
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source  is  bi-directional.  iDSIM2  also  detects  subcircuits  which  have  the  same  circuit 
configuration  and  builds  a  two-level  hierarchy.  Although  iDSIM2  also  includes  other  higher 
level  simulation  options  such  as  switch-level  and  logic-level  simulation,  only  circuit  and  fast¬ 
timing  are  involved  in  reliability  analysis  because  of  the  accuracy  requirements. 

The  use  of  the  device  degradation  models  must  be  restricted  to  shorter  time  intervals  in 
which  the  change  of  stress  conditions  may  be  neglected.  It  is  assumed  that  the  degradation 
model  can  be  applied  to  predict  the  hot-carrier  induced  device  degradation  accurately  within  a 
time  interval  of  length  T,.  The  goal  is  to  predict  the  amount  of  device  and  circuit  perfor¬ 
mance  degradation  for  the  operation  time  interval  T2  >  T,.  The  period,  Tq,  of  the  input 
waveform  is  assumed  to  be  Tq  «c  Ti,T2. 

Simulation  procedure  ; 

(i)  Simulate  the  circuit  for  one  period  (To)  and  determine-  the  stress  conditions  (terminal  vol¬ 
tage  waveforms)  associated  with  each  nMOS  transistor. 

(ii)  Using  the  hot-carrier  degradation  models,  determine  the  amount  of  degradation  that  each 
transistor  will  experience  at  the  end  of  the  time  interval  T]. 

(iii)  Update  the  damage  parameters  of  each  transistor  to  specify  the  hot-carrier  degradation 
sustained  at  the  end  of  T^. 

(iv)  Return  to  (i)  and  repeat  the  procedure  k  times  such  that  T2  =  k  T). 

It  can  be  seen  that  the  results  of  each  circuit  simulation  cycle  (for  one  period  Tq)  are 
used  to  estimate  the  amount  of  hot-carrier  induced  damage  each  transistor  will  experience  dur¬ 
ing  a  time  interval  of  length  T,. 

After  each  one-cycle  simulation,  the  damage  level  of  each  transistor  is  extrapolated  to  the 
previously  specified  time  period,  and  the  parameters  are  updated  accordingly.  Based  on  this 
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information,  the  user  may  want  to  switch  those  subcircuits  which  have  no  degradation  to  fast- 
dming  in  order  to  increase  speed.  The  user  can  also  change  the  operation  duration  if  the  pre¬ 
vious  one  is  too  long  or  too  short,  and  re-run  the  simulation.  With  this  mixed-mode  simula¬ 
tion  capability,  iDSIM2  can  handle  large  circuits  with  computational  efficiency. 

The  use  of  macromodels  is  also  being  considered  for  representing  damaged  subcircuits  in 
circuit-level  simulation  using  iDSIM2.  This,  combined  with  the  use  of  fast-timing  macromo¬ 
dels  for  undamaged  subcircuits,  will  further  improve  the  computational  efficiency  while 
preserving  the  required  accuracy. 

A  824-transistor  combinational  circuit  has  been  simulated  to  demonstrate  the  large-scale 
simulation  capability  of  iDSIM2.  The  simulation  cycle  length  in  this  example  is  Tq  =  84  ns. 
Following  the  one-cycle  simulation,  hot-carrier  induced  degradation  of  the  nMOS  transistors 
are  estimated  for  time  interval  of  Tj  =  5x10^  seconds.  After  updating  the  damage  parameters, 
the  simulation  is  repeated  as  described  earlier.  The  total  CPU  time  for  performing  this  simu¬ 
lation  on  a  Sun  3/160  is  48  minutes.  The  distribution  of  the  simulated  hot-carrier  damage  lev¬ 
els  in  nMOS  transistors  is  given  in  Fig.  14.  The  second  peak  of  the  distribution  graph  indi¬ 
cates  the  transistors  with  extensive  hot-carrier  induced  damage,  which  can  influence  the  circuit 
performance  characteristics.  Thus,  any  gate  containing  at  least  one  extensively  damaged 
transistor  must  be  classified  as  a  "damaged  gate." 

V.  Design  of  Reliability  Tester 

The  reliability  simulation  tools  developed  in  this  task  need  to  be  verified  and  calibrated 
using  experimental  data  on  hot-carrier  degradation  of  MOS  circuits  under  realistic  operating 
conditions.  To  achieve  this  goal,  a  tester  chip  containing  a  42-stage  inverter  chain  ,ias  been 
designed  at  the  University  of  Illinois,  and  fabricated  by  Sematech  in  the  1990-1991  period.  In 
the  tester  chip,  some  transistors  have  been  designed  to  be  individually  controllable  and 
observable,  for  accurate  determination  of  device  characteristics  before  and  after  hot-carrier 
stress.  The  circuit  can  be  stressed  under  different  dynamic  operating  conditions  for  different 
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lengths  of  operation.  The  changes  in  the  transistor  characteristics  can  be  monitored  at  certain 
intervals  during  long-term  stressing  in  order  to  determine  the  amount  of  hot-carrier  stress 
experienced  by  each  device. 

The  transistors  to  be  accessed  externally  are  organized  into  two  separate  test  ceils,  each 
of  which  contains  two  CMOS  inverters  in  series.  The  first  test  cell  constitutes  stage  15  and 
stage  16  of  the  inverter  chain,  whereas  the  second  test  cell  constitutes  stage  31  and  stage  32. 
Except  for  the  test  cells,  the  (W/L)-ratios  of  nMOS  and  pMOS  transistors  are  (14.4/1.2)  and 
(29.4/1.2),  respectively.  The  test  cells  contain  one  inverter  with  the  same  (W/L)-ratios  as  the 
regular  stages,  and  one  inverter  with  (W/L)-ratios  of  (9.6/1. 2)  and  (14.4/1.2).  Two  versions  of 
the  test  chip  have  been  manufactured  using  the  Phase-2  and  Phase-3  design  rules  developed 
by  Sematech. 

As  indicated  on  the  test-cell  circuit  diagram  in  Fig.  15,  8  pads  are  dedicated  to  each  test 
cell  so  that  ail  four  transistors  in  a  cell  are  independently  accessible.  This  feature  allows  the 
accurate  measurement  of  device  model  parameters  of  each  transistor  under  various  bias  condi¬ 
tions,  before  and  after  hot-carrier  stress.  The  transistors  in  the  test  cells  can  be  stressed  under 
dynamic  circuit  operating  conditions,  or  under  externally  applied  DC-bias  conditions  to 
accelerate  the  degradation  of  individual  devices.  In  addition  to  the  3  pads  for  Vss  and  2  pads 
for  Vdd,  one  pad  each  is  provided  for  signal  input,  41st  stage  output  and  42nd  stage  output. 
The  41st  stage  output  can  be  used  to  convert  the  inverter  chain  circuit  into  a  ring  oscillator  by 
an  external  connection. 

VI.  Summary 

During  the  last  project  period,  we  have  developed  new  simulation  capabilities  for  relia¬ 
bility  analysis  of  MOS  circuits  with  emphasis  on  hot-carrier  damages.  In  particular,  both 
iSMILE  and  iDSIM2  have  been  used  to  develop  and  implement  new  reliability  models  of 
nMOS  transistors. 
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The  physical  one-dimensional  model  for  hot-carrier  damaged  nMOS  transistor  has  been 
shown  to  explain  the  unique  asymmetry  in  I-V  characteristics  due  to  local  damages  in  the 
drain  of  the  channel,  besides  the  changes  in  the  transconductance  and  the  mobility  degrada¬ 
tion.  This  model  can  be  incorporated  in  SPICE-like  circuit  simulators  as  demonstrated 
through  iSMILE  which  is  a  model-independent  circuit  simulator  that  enables  rapid  model 
development  Simulation  results  have  shown  strong  agreements  with  measured  data  on 
transistor  characteristics,  which  points  to  potential  matches  at  circuit  level.  However,  the  cir¬ 
cuit  level  simulation  capability  remains  to  be  verified  in  the  future  due  to  insufficient  data. 

To  extend  the  circuit  simulation  capability  to  large  MOS  circuits,  we  have  also 
developed  a  two-tier  mixed-level  reliability  simulation  approach  using  the  iDSIM2  program. 
Using  iDSIM2,  we  have  been  able  to  simulate  a  824-transistor  combinational  circuit  on  a 
Sun- 3/1 60  workstation.  This  is  a  major  speed  improvement  over  iSMILE-based  simulation 
and  allows  in-depth  reliability  analysis  of  large  MOS  circuits. 

In  an  effort  to  verify  the  accuracy  of  circuit  simulation  using  iSMILE  or  iDSIM2,  we 
have  designed  a  CMOS  tester  wherein  I-V  characteristics  of  some  transistors  can  be  measured 
individually  before  and  after  hot-carrier  stress.  This  tester  is  being  fabricated  by  Sematech. 
Future  experiments  with  the  fabricated  tester  chips  should  be  helpful  in  improving  the  model 
and  the  simulation  capability. 

Besides  the  planned  experimental  verification  tasks,  other  basic  tasks  need  to  be  investi¬ 
gated.  To  mention  a  few,  development  of  models  for  pMOS  transistors  and  lightly  doped 
drain  (LDD)  structures  would  be  important  in  view  of  technology  trends.  Furthermore,  to 
cope  with  the  ever-increasing  level  of  device  integration,  a  faster  simulator  that  can  handle 
several  hundreds  of  thousands  of  transistors  needs  to  be  developed. 
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PARAMETER 


VALUE 


W  (cm) 

10.0x10-^ 

L  (cm) 

i.oxicr* 

Ho  (cm^/y^s) 

678.5 

(V) 

-0.8385 

<J>P  (V) 

0.354 

(cm~^) 

1.809x10^^ 

N„  (cm’^) 

1.0x10^° 

tax 

224x10'* 

a 

0.3188 

P 

0.4982 

T1 

0.0802 

Ec  (Vicm) 

2.36x10'^ 

Table  1 
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1.18 


POSITION  ALONG  CHANNEL  0*fn) 


Fig.  1.  Measured  interface  charge  distribution  profile  near  drain  [12]. 
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Fig.  2.  Oxide  charge  distribution  profile  used  in  model  derivation. 
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n  Damaged  oxide  region 

^  Gaussian  surface  enclosing  mobile  and  bulk  charges 
Channel  region  where  GCA  is  valid 


Fig.  3.  Pseudo  two-dimensional  nxxieling  of  the  damaged  MOSFET. 
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DRAIN  CURRENT  (A) 


Fig.  4.  Comparison  of  measured  and  simulated  Id-Vq  characteristics  of  the  nMOS  transistor 
before  and  after  hot-carrier  stress.  Stress  conditions  ;  Vq  =  2.6  V,  Vp  =  6.5  V, 
*  423  h. 

O  :  Measured  drain  current  before  stress 
A  :  Measured  drain  current  (tfter  stress 
Solid  lines  ;  Simulated  drain  current  before  stress 
[>otted  lines  :  Simulated  drain  current  after  stress 


32 


DRAIN  CURRENT  (mA) 


1.4 


DRAIN  VOLTAGE  (V) 


Fig.  7.  Forward  and  reverse  Id-Vq  characteristics  of  the  nMOS  transistor  before  and  after 
hot-carrier  stress.  Stress  conditions  :  Vc  =  3  V,  Vp  =  8  V,  t„re„  =  5  x  lO'*  s.  Channel 
length  L  =  1.7  nm  (7], 

O  ;  Measured  (forward  and  reverse)  drain  current  before  stress 
A  :  Measured  forward  drain  current  after  stress 
□  :  Measured  reverse  drain  current  after  stress 
Solid  lines  :  Simulated  drain  current  before  stress 
Dashed  lines  ;  Simulated  drain  current  after  stress 
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0.00015 


GATE  VOLTAGE  (V) 


Fig.  8.  Simulated  Iq-Vg  characteristics  of  the  nMOS  transistor  with  different  levels  of  hot- 
cairier  damage  N;,. 

Solid  line  ;  Undamaged  device  (N;,  *  0). 

Dashed  line  :  Nj,  =  0.5x10*^  cm~^ 

Dash-dotted  line  :  N,,  =  l.OxiO*^  cm“^. 

Dotted  line  :  Nj,  =  l  .Sxlo'^  cm“^. 
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TIME  (ns) 


Fig.  9 


Simulated  operation  of  the  7-stage  ring  oscillator  before  (solid  lines)  and  after  (dotted 
lines)  hot-carrier  induced  oxide  damage. 

(a)  Locally  damaged  nMOS  transistors 

(b)  Uniformly  damaged  nMOS  tranristors 


VOLTAGE  (V) 


Fig.  10.  Simulated  operation  of  the  cross-coupled  latch  circuit  (inset),  with  locally  damaged 
transistors  (solid  lines)  and  with  uniformly  damaged  transistors  (dotted  lines). 
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Fig.  11.  nMOS  Pass- transistor  circuit. 
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VOLTAGE  (V) 


TIME  (s) 


Fig.  12.  Simulated  input  and  output  voltage  waveforms  of  the  nMOS  pass-transistor  circuit. 

Solid  line  :  Before  degradation. 

Dashed  line  :  After  degradation. 
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Fig.  13.  Structure  of  the  reliability  simulation  procedure. 
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0.7  0.8  0.9  1  2  3  4  5  6  7  8  9  10  20  30 

Damage  level  (Nit)  (x  lrf®/cm^ 

Fig.  14.  Distribution  of  simulated  hot-carrier  damage  levels  in  nMOS  transistors  after  SxlO^ 
seconds  of  dynamic  operation. 
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MISSION 

OF 

ROME  LABORATORY 

Rome  Laboratory  plans  and  executes  an  interdisciplinary  program  in  re¬ 
search,  development,  test,  and  technology  transition  in  support  of  Air 
Force  Command,  Control,  Communications  and  Intelligence  (C  I)  activities 
for  all  Air  Force  platforms.  It  also  executes  selected  acquisition  programs 
in  several  areas  of  expertise.  Technical  and  engineering  support  within 
areas  of  competence  is  provided  to  ESD  Program  Offices  (POs)  and  other 

O 

ESD  elements  to  perform  effective  acquisition  of  C  l  systems,  bi  addition, 
Rome  Laboratorys  technology  supports  other  AFSC  Product  Divisions,  the 
Air  Force  user  community,  and  other  DOD  and  non-DOD  agencies.  Rome 
Laboratory  maintains  technical  competence  and  research  programs  in  areas 
including,  but  not  limited  to,  communications,  command  and  control,  battle 
management,  intelligence  information  fsrocessing,  computational  sciences 
and  software  producibility,  wide  area  surveillance/sensors,  signal  proces¬ 
sing,  solid  state  sciences,  photonics,  electromagnetic  technology,  super¬ 
conductivity,  and  electronic  reliability/maintainabUity  and  testability. 


